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ABSTRACT 

We present X-ray spectral analysis of five Chandra and XMM-Newton obser- 
vations of the gravitationally-lensed blazar PKS 1830—211 from 2000 to 2004. 
We show that the X-ray absorption toward PKS 1830—211 is variable, and the 
variable absorption is most likely to be intrinsic with amplitudes of ~ 2 x 10^^- 
30 X lO^^cm"^ depending on whether or not the absorber is partially covering 
the X-ray source. Our results confirm the variable absorption observed previ- 
ously, although interpreted differently, in a sequence of ASCA observations. This 
large variation in the absorption column density can be interpreted as outfiows 
from the central engine in the polar direction, consistent with recent numerical 
models of infiow/outfiows in AGNs. In addition, it could possibly be caused by 
the interaction between the blazar jet and its environment, or the variation from 
the geometric configuration of the jet. While the spectra can also be fitted with 
a variable absorption at the lens redshift, we show that this model is unlikely. 
We also rule out the simple microlensing interpretation of variability which was 
previously suggested. 

Subject headings: 



1. Introduction 



PKS 1830—211 (jPramesh Rao fc Subrahmanvan lll988l : ISubrahmanyan et al. lll990l : lJauncey et al. 
1991 ) consists of two Zg = 2.507 (Lidman et al. I 1999 ) blazar images se parated by 1" 



and lensed by a = 0.886 ( Wiklind fc Combes 19961 : Gerin et al. 1997 ) spiral galaxy 
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1988 



Lovell et al 



( Winn et al. 112002). PK S 1830—211 w as observed in the radio (e.g., Pramesh Rao fc Subrahnaanvan 



19981) ■ infrared (e.g 



20021), X-rays (IMathur fc Nair 



1997 



Oshima et al 



2001 



Courbin et al. 


2002; 


Winn et al. 


de Rosa et al. 


2005: 


Dai et al. 



2006f). and Ganima-ra ys (IMattox et al.lll997l ). and the spectral energy distribution of PKS 1830—211 
(jde Rosa et al. II2005I ) shows two emission bumps with one in the infrared and the other be- 
tween the hard X-ray and Gamma-ray bands. 

The gravitational lens PKS 1830—211 is complicated in many aspects, a nd one of them is 



the X -ray absorption. The X-ray a bsorption was first d etected by ROSAT (IMathur fc Nair 



19971 ). Later ASCA observations (jOshima et al.ll200ll ) show that the X-ray absorption is 
variable; these authors favored a model in which microlensing was the cause of variability. 
Recently, PKS 1830—211 was observed by Chandra and XMM-Newton for five epochs, which 
enabled us to study the nature of the X-ray absorption in detail combining the angular 
resolution of Chandra and the large collecting area of XMM-Newton. Since the two images 
are resolved by Chandra., we can test the microle nsing mode l predi ctions and study the X-ray 
absorption separately for the two lensed images. iDai et al. I (120061 ) found that the differential 



absorption at the lens galaxy between the two images is AiYns^A = 1.8lQ;gX 10^^ cm ^. In this 



paper, we study the time evolution of the absorption in the system from the five Chandra and 
XMM-Newton observations and confirm that the X-ray absorption toward PKS 1830—211 is 
variable. Moreover, we rule out the microlensing interpretation of the observed variability. 



2. Observations and Data Reduction 

We observed PKS 1830—211 twice with Chandra and three times with XMM-Newton 
from 2000 to 2004. The details of the observations are listed in Table [U The Chandra 
data were reduced with the CIA03 . 2 software tools provided by the Chandra X-ray Center 
(CXC) following the standard threads on the CXC websiteE Only events with standard 
ASCA grades of 0, 2, 3, 4, and 6 were used in the analysis. We improved the image quality 
of the data by removing the pixel randomization applied to the event p ositions by the 



standard pipeline. In a ddition, we applied a sub-pixel resolution technique (ITsunemi et al. 



200ll ; iMori et al. 1 120011 ) to the events on the S3 chip of ACIS where the quasar images are 
located. This allowed us to resolve the two lensed images of the blazar. The XMM-Newton 
data were reduced using the standard analysis software SAS6 . 0. We used the tasks epchain 
and emchain from SAS to reduce the PN and MOS data, and filtered the events with patterns 
< 4 and < 12 for the PN and MOS chips, respectively. 



^The CXC website is at http://cxc.harvard.edu/ 
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3. Spectral Analysis 



We fitted the spectra of PKS 1830-211 using XSPEC Vll.3.1 flArnaudl Il996[ ) in the 
0.35-8 keV observed energy range for Chandra spectra and in the 0.35-10 keV range for 
XMM-Newton sp ectra. In all of our mode ls, we fixed the Galactic absorption at iVn = 
0.22 X 10^2 cm-2 Jpickev fc Lockman lEoQol ). 



3.1. Chandra Spectra of Individual Images 

We analyzed the Chandra spectra of individual images A and B. We modeled the spectra 
with a power law modified by the Galactic absorption and the absorption at the lens redshift; 
the fitting results are listed in Table [2J We experimented with two models where the power 
law photon indices for the two images were allowed to be different (Model 1) and where they 
were constrained to be the same (Model 2). We assumed that the excess absorption (above 
Galactic) arises in the lens galaxy at z=0.886 (we test this assumption later). In both the 
models, the column densities at the two Chandra epochs for image B are consistently larger 
than those of image A. In addition, the absorption in the first epoch is consistently larger 
than that in the second epoch for both images, except for image B in Model 1 where it is 
similar within errors. The photon indices obtained from Model 1 fits are similar for the two 
images and in both epochs, partially due to the large error bars on the parameter caused 
by the low signal-to-noise ratio of the individual spectra. The photon indices obtained from 
Model 2 show ~ la difference from epoch to epoch, which again is not significant. The 
difference between the absorption in the two images occurs simply because t he lines of sight 



of the two images intersect different parts of the lens galaxy (jPai et al. II2006I ). The difference 
between the absorption at the two epochs indicates that the X-ray absorption is variable. 
The origin of the absorption variability could be either at the lens redshift or at the source 
redshift. We note that the differential absorption between images B and A is similar in the 
two epochs, especially in Model 2. This is suggestive of absorption variability occurring at 
the source redshift. 



3.2. Chandra and XMM-Newton Spectra of Combined Images 

We then analyzed the spectra of the combined images AB for the Chandra and XMM- 
Newton observations. We first fitted the spectra of PKS 1830—211 of the five epochs with a 
power-law model modified by neutral absorption from the Milky Way and the lens galaxy. 
The fitting results are listed in Table [3] (Model 3). We note that the absorption at the 
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lens redshift for this model should be treated as an averaged absorption for the two lines of 
sight. We also fitted the co-added three XMM- Newton spectra to obtain a higher signal-to- 
noise ratio spectrum, and the results are also listed in Table [31 The higher S/N of XMM- 
Newton spectra allows better constraints on the power-law spectral index. We find that 
the spectral index varies between the Chandra and the XMM-Newton observations, though 
variations within the Chandra epochs and XMM-Newton epochs are smaller. Comparing the 
fitting results of Model 3 and Model 2, the spectral index obtained by fitting the combined 
images AB for the Chandra observations are consistent with results from individual spectral 
fits in § 13. 1[ The absorption at the lens is also variable, especially when comparing the 
second Chandra observation with other epochs. We further test the variability of the spectra 
by comparing a model with no spectral variability (with the exception of normalization) 
and another model with variable absorption and spectral index. We found that the model 
with the spectral variability produced a better fit with a Ax^ improvement of 56.1 (a null 
probability of 10~^ by the F-test) when jointly fitting the five X-ray observations. When 
we only consider the three XMM-Newton observations, the corresponding improvement is 
14.5 with a null probability of 0.0097 given by the F-test. The largest absorption variation is 
AA^H = (0.7±0.3) X 10^^ cm~^ between the second Chandra epoch and t he first XMM-Newto n 



epoch. Although the variable spectral index in blazars is common (e.g.. lFoschini et al.ll2006l ) 



the absorption at the lens galaxy is unlikely to vary on the time scales of years. As discussed 
in § 13. H it is more likely that the variable absorption component is at the source redshift 
close to the AGN, where a short time-scale variability is possible and may be expected. 

The next model we tried is a power-law modified by three absorption components, the 
Galactic absorption and the absorptions at the lens and the source. We fixed the absorption 
at the lens as Ah = 1.7 x 10^^ cm~^, the smallest absorption column density detected in 
the five epochs from the previous model. This is the largest absorption column density at 
lens we can set to ensure no absorption variability at the lens galaxy. The fitting results 
are listed in Table [3] (Model 4). Again we detected variability of both the photon index and 
the absorption at the source redshift. Using the same test that we described in the previous 
model, we found that the model allowing spectral variability improved the Ax^ by 59.8 with a 
null probability of 5 x 10"^ for jointly fitting the five X-ray observations, and an improvement 
of 14.7 with a null probability of 0.01 for jointly fitting the three XMM-Newton observations. 
The largest absorption variation for this model is A Ah = (3.5 ± 0.7) x 10^^ cm~^ between 
the second Chandra epoch and the first XMM-Newton epoch. The variable absorption in 
this model can be naturally associated with outflows from the central engine. Although, 
this model produces comparable fits for the two Chandra spectra compared with Model 3, 
where there is no absorption at the source redshift, the fits for the XMM-Newton spectra 
with a higher S/N are worse than Model 3 with Ax^ increases of 5.8, 37.2, 13.0 for the 
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three XMM-Newton spectra. We also tested whether fitting the absorption component at 
the source with a warm absorber or a partially covering absorber would improve the fit. 
While the warm absorber model produces a worse fit, the partial covering model (Table [31 
Model 5) produces comparable fits to Model 3. For the partial covering model, the model 
allowing spectral variability produces improvement of the Ax^ by 79.7 with a null probability 
of 1 X 10~^^ for jointly fitting the five X-ray observations, and an improvement of 14.1 with a 
null probability of 0.01 for jointly fitting the three XMM-Newton observations. The largest 
absorption variation for this model is AA^h = (29 ±7) x 10^^ cm~^, again, between the second 
Chandra epoch and the first XMM-Newton epoch. 

As a final improvement to our model, we assume that the lines of sight of the two lensed 
images pass through different locations in the lens galaxy and thus have different absorption 
column densities. This differential absorption is measured as AA^hb.a = l-8-0 6 ^ ^^^^ cm~^ 



( iDai et al. II2006I ). Therefore, we have a model composed of two power-law components and 
each of them is modified by three absorptions from Galactic, lens, and source. The two 
power-law components representing the two lensed images have the same photon index but 
different normalizations. The normalization ratio between images A and B is constrained 
as i? = 1.03 by taking the hard X-ray flux ratios (3-8 keV) from the Chandra observation. 
The Galactic absorption and the absorption at the source are the same for the two power- 
law components, and the absorbers at the source are assumed to be partially covering the 



continuum. The absorber column densities at t 



components by the amount given by lDai et al. I (120061 ). We flt the two Chandra spectra and 



le len s are different for the two power-law 



the co-added XMM-Newton spectrum simultaneously with this model. We constrained that 
the absorption at lens did not vary. The fitting results are listed in Table [3] (Model 6), and 
the spectra are shown in Figure [H The fitting results indicate that the absorption at the 
source has varied by roughly (39 ± 7) x lO^^cm"^. 



4. Discussion 

We detected the variable absorption, variable power-law photon index, and variable flux 
for the blazar PKS 1830—211. The observed variability between Chandra and XMM-Newton 
spectra could be, in part, due to imperfect cross-calibration between the two. However for 
high S/N spectra, the cross-calibration between Chandra and XMM-Newton can yield a 
spectral difference of AF ~ 0.03|^, whereas our measured differences are between AF ~ 
0.15-0.3. In addition, the higher S/N XMM-Newton spectra show a statistically signiflcant 



^http: / /xmm.esac.esa.int/external/xmni_sw_cal/calib/cross_cal/index.php 
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variability within the three XMM-Newton observations (Ax^ improvement of 14.5, 14.7, and 
14.1 for Models 3, 4, and 5 with spectral variability). We also found the spectral variability 
between the two Chandra observations, where the second Chandra epoch has consistently 
lower absorption column densities than the first observation. This is better demonstrated 
in § 3.1, where we analyzed the individual spectra of each image. We plot the variations of 
absorption in Figure [21 photon index variations in Figure [3l and the photon index against 
unabsorbed flux in Figure |H It is well known that blazars are variable, so the observed 
spectral variability is not surprising in itself. We also found that the photon index and flux 
are correlated (Figure H]). This relation has been observed in several Gamma-ray loud AGNs 
(Foschini et al. 2006). PKS 1830—211 provides an additional example of a source following 
this correlation. 

The variable absorption is more intriguing. T he variable absorptio n toward PKS 1830—211 



was detecte d previously w i th AS CA observations (lOshima et al.ll200ll ). although interpreted 



differently. lOshima et al.l (120011 ) obtained a variation amplitude of ~ 0.5 x lO^^cm"^ when 
modeling the absorption at the lens, which is consistent wit h our measurernent o f ~ 0.6 x 
lO^^cm"^ (Model 3). We note that the variation detected by lOshima et al.l (1200 ll ) is on the 



timescale of 10 days while that found by us is on the timescale of years. It is possible that 
the variable absorption is at the source redshift. Depending on the nature of the absorber, 
we obtained variations of ~ 2 x lO^^cm"^ and ~ 30 x lO^^cm"^ for a fully covering absorber 
(Model 4) and a partially covering absorber (Model 5, 6), respectively. The other possibilities 
are microlensing models. 



Oshima et al.l (120011 ) proposed a microlensing interpretation to explain the apparent 
variable absorption. The basic idea is that the absorption does not vary but is different 
for the two images, and the microlensing produces flux changes between the two images 
which appears as variable absorption. As we discussed in the spectral fits to individual 
Chandra images (§ 13.11) . the absorption varies in individual images between the two epochs, 
which cannot be attributed to this simple microlensing model. For this reason, we rule 
out the simple microlensing model as the cause of variable absorption in PKS 1830—211. 
More complicated microlensing models are required which must involve significant emission 
size differences between the soft and hard X-ray regions to interpret the data. However, 
even these models will have difficulty explaining the similar differential absorptions between 
images A and B for the two Chandra epochs. Microlensing will only produce uncorrelated 
changes between the images. 

Variable absorption thus appears to be a robust result, independent of model. The vari- 
ability can either come from the lens or from the source redshift. If the variable absorption 
is at the lens, a variability amplitude of ~ 0.6 x lO^^cm"^ can sufficiently parameterize the 
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Chandra and XMM-Newton spectra without invoking any absorption at the source. Re- 
cently, variable Mg I I absorption has been detected in an intervening absorber in front of 
a Gamma-ray burst (jHao et al. 2006). The obs erved Mg II column density changes as the 
burst expands geometrically (IFrank et al. II2006I ). Similarly, if the size of the blazar changes 
with luminosity, it may produce variable absorption at the lens. Another mech anism that 
can produce variable absorption at the lens galaxy is proposed by iDong I (120061 ) . where mi- 
crolensed images, slightly shifted from the original image position, intersect slightly different 
positions at the lens, producing variable absorption. However, in PKS 1830—211, the differ- 
ential absorption of the two images in the two Chandra epochs is similar, arguing against 
the variability at the lens. It is far more likely that the variable absorption is intrinsic to the 
blazar, simply because there are other AGNs with variable intrinsic absorption, including 
blazars (see below). 

While the absorption toward the gravitational lens PKS 1830—211 is very complicated, 
with the sequence of Chandra and XMM-Newton observations, a conver ging picture is emerg- 



mg. 



Besides the Galactic absorption of A'^h = 22.27 x 10 cm ( iDickey fc Lockman 



19901), there is also an intrinsic absorption and absorption from the lens galaxy. The 
absorption column densities in the lens galaxy are different for the two lensed images by 



^NuB,A = l-8-a6 ^ 10 cm ^ ( iDai et al. II2006I ). and this is consistent with the differential 



extinction measurement between the two images and a Galactic dust-to-gas ratio. The am- 
plitude of the variabihty of the intrinsic absorption is model dependent with amplitudes of 
~ 2 X 10^2-30 X 1022cm-2. 

Recently, the int rinsic absorption for blazars has been rep orted in several cases , such 
as GB B1428+4217 JWorslev et al.ll2004ah . PMN J0525-3343 JWorslev et al.ll2004bh . and 
RBS 315 (Piconcelli & Guainazzi 2005; see also Tavecchio et al. 2007). In this paper, we 
present variable intrinsic absorption toward gravitationally-lensed blazar PKS 1830—211. 
The large variation in the absorption column density can be interpreted as outflows similar 
to those detected in BAL QSOs. However, as PKS 1830—211 is a blazar, this outflo w must 



be in the polar direction. Recent numerical simulations of accretion flows in AGNs (jProga 



20061 ) discussed the cases of polar outflows; our observations of PKS 1830—211 might be 
providing an example in support of such models. In addition, it is also possible that the 
intrinsic absorption variation is caused by the interaction between the blazar jet and its 
environment, or the variation from the geometric configuration of the jet. More intense 
monitoring of this system or further variability studies using a large sample of blazars are 
needed to better constrain the nature of their X-ray spectral variation. 
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Table 1. List of X-Ray Observations of PKS 1830-211 



Date 


Telescope 


Grating? 


Exposure Time (sec) 


2000-06-26 


Chandra 


HETGS 


47471 


2001-06-25 


Chandra 


HETGS 


51219 


2004-03-10 


XMM-Newton 


none 


8328 


2004-03-24 


XMM-Newton 


none 


27004 


2004-05-05 


XMM-Newton 


none 


21356 



Table 2. Spectral Analysis for Individual Chandra Spectra of PKS 1830—211. 



Date Image ArH(z=0.886), 10^2 cm'^ r X^idof) 



Model 1: An absorbed power-law 



2000-06-26 


A 


1.6 ±0.4 


1.09 ± 0.09 


100.5(108) 


2000-06-26 


B 


3.5 ± 1.1 


1.13 ± 0.15 


129.6(106) 


2001-06-25 


A 


0.8 ±0.5 


0.96 ± 0.09 


140.7(120) 


2001-06-25 


B 


3.3 ± 1.0 


1.11 ± 0.12 


113.8(116) 




Model 2; 


An absorbed power- 


-law with = ^B- 




2000-06-26 


A 


1.6 ±0.4 


1.10 ±0.07 


229.6(215) 


2000-06-26 


B 


3.4 ±0.7 






2001-06-25 


A 


1.0 ±0.4 


1.02 ± 0.06 


255.2(237) 


2001-06-25 


B 


2.8 ±0.6 







Note. — Galactic ATh is fixed at ATr = 22.27 x 10^° cm-^ (Dickey & Lockman 
1990) for all models. 



Table 3. Spectral Analysis of Combined Images of PKS 1830—211. 



Chandra XMM-Newton 



Parameters 


I 


II 


I 


II 


III 


XMM-Co-Added 


Date 


2000-06-26 


2001-06-25 


2004-03-10 


2004-03-24 


2004-05-05 




Model 3: A power-law with absorption at the lens and the Galaxy 


Nji{z = 0.89,1022 cm-2) 

r 

Unabsorbed Flux ^ 
xHdof) 


2.27 ±0.26 

1.07 ±0.05 
0.92 ±0.02 
365.6(305) 


1.70 ±0.22 

0.99 ±0.04 
0.97 ± 0.03 
365.6(337) 


2.42 ±0.12 

1.18 ±0.02 
1.50 ±0.03 
151.8(152) 


2.37 ±0.07 

1.14 ±0.01 
1.42 ± 0.02 
447.3(396) 


2.29 ± 0.07 

1.17 ± 0.01 
1.28 ±0.02 
314.3(301) 


2.34 ±0.04 

1.16 ± 0.01 
1.38 ± 0.01 
843.7(743) 


Aid 


(1(4 i: Same as 


Ai(>(l(-1 :l Willi 


addilioiial ai>^ 


aii])l i(>].i at liie 


source'' 




Nb_{z = 2.51,1022 cm-2) 

r 

Unabsorbed Flux 
xHdof) 


1 

-'-■^^-1.19 

1.02 ± 0.05 
0.91 ± 0.04 
369.9(305) 


< 0.38 
0.99 ±0.03 
0.97 ±0.03 
365.6(337) 


3.47 ±0.60 
1.17 ±0.02 
1.50 ±0.03 
157.6(152) 


2.80 ± 0.33 
1.12 ±0.01 
1.41 ± 0.02 
484.5(396) 


2.68 ±0.36 
1.16 ±0.01 
1.28 ±0.02 
327.3(301) 


2.80 ± 0.22 
1.14 ±0.01 
1.38 ±0.01 
899.3(743) 


Model 5: 


Same as Model 4, but with a 


partially covering absorber at the source'' 




Nu{z = 2.51, 1022 cm-2) 
Covering Factor 

r 

Unabsorbed Flux 


< 4.6 
0.44 (iixcd) 
1.00 ±0.04 
0.91 ±0.04 
370.6(305) 


< 0.7 
0.44 (fixed) 
0.98 ±0.03 
0.97 ± 0.04 
365.9(337) 


29. oil;" 
0.50 ± 0.05 
1.32 ±0.05 
1.70 ± 0.06 
137.4(151) 


0.43 ± 0.03 
1.23 ± 0.03 
1.53 ±0.04 
447.4(395) 


19 5+'^-l 
0.40 ± 0.04 
1.23 ±0.03 
1.36 ± 0.04 
315.6(300) 


25.0 ± 2.5 
0.44 ± 0.02 
1.26 ± 0.02 
1.51 ±0.03 
827.0(742) 



Model 6: Same as Model 5, but with two diflferent power-law components'' 



NuA = 0-89, 10^2 crn-2) I.57 ± 0.07 

ArH(2 = 2.51, 1022 cm-2) < Q 7 ... ... ... 39 ±7 

Covering Factor 0.23 ± 0.04 

r 1.07 ±0.03 1.07 ±0.02 ••• ■■■ ••• 1.27 ±0.03 

X^dof) 1488.8(1367) 



Note. — Galactic Nh is fixed at A^h = 22.27 x 10^° cm-2 pj^key k Lockman 1990) for all models. 
''The unabsorbed fluxes in all models are calculated between 0.4-8 keV and in units of 10-'2erg cm-2 g-i 
''The absorption at lens is fixed at Nn = 1.7 x IO22 cm-2 j^j. ^j^jg model. 

''The absorption of the lens galaxy for image B is constrained as Nug = + 1-8 x IO22 cm-2. 



- 12 - 




Observed — Frame Energy (keV) 



Fig. 1. — XMM-Newton and Chandra spectra of PKS 1830—211. The top spectrum is a 
combined spectrum from three XMM-Newton observations, and the bottom two spectrum 
is from the two Chandra observations. The spectra are fitted simuhaneously with absorbed 
power-law models as described in Model 6 of Table [31 
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Fig. 2. — The absorption column densities at the redshift of the source {z — 2.51) of 
PKS 1830—211 measured from Chandra and XMM-Newton spectra assuming a partial cov- 
ering model (filled circles) and a fully covering absorption model (squares). 
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Fig. 3. — The photon indices measured for PKS 1830—211 from Chandra and XMM-Newton 
observations. 
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Fig. 4. — The photon indices versus flux for PKS 1830—211 from Chandra and XMM-Newton 
observations. 



